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Abstract 

The oxidation of SiC-protected carbon/carbon through machined slots and naturally occurring craze cracks in the 
SiC was studied. The slot and crack geometries were characterized, and the subsurface oxidation of the 
carbon/carbon substrate at temperatures of 1000 to 1300 °C in air was assessed using weight change, x-ray 
computed tomography, and optical microscopy of sections. Rate constants were derived from these measurements 
and compared with a two-step diffusion control model of carbon oxidation. Oxidation kinetic measurements on both 
the specimens with machined slots and with naturally occurring craze cracks showed good agreement with the 
model. 
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Introduction 

Reinforced Carbon/Carbon 

Fig. 1 is a schematic of reinforced carbon/carbon (RCC) with a SiC conversion coating used on the Space Shuttle 
Orbiter nose cap and wing leading edges. The material is made with a two-dimensional layup of carbon-carbon 
fabric with repeated applications of a liquid carbon precursor to fill voids [1]. The resultant high strength and light 
weight make it an ideal aerospace material; however, oxidation is a major concern. An oxidation protection system 
is based on a SiC conversion coating. Because of the difference in coefficient of thermal expansion (CTE) of the SiC 
coating and carbon/carbon substrate, the SiC coating shrinks more than the underlying carbon/carbon on cooldown 
from the coating application temperature. This leads to vertical cracks in the coating, and these cracks are pathways 
for oxygen to reach the carbon/carbon substrate. Actual RCC used on the orbiters contains glass sealants to plug the 
cracks; however, this work examines carbon/carbon with only SiC. 

The crack morphologies in a similar system have been thoroughly characterized in an elegant study by 
Yurgartis et al. [2]. Both surface (plan) views and section views of the material were taken to obtain information 
such as crack spacing, crack width, and crack morphology. Such parameters were obtained for RCC in this study 
and prove critical to modeling oxidation. Buchanan and Little [3] have also characterized the crack pattern in a range 
of SiC coatings on carbon/carbon at various deposition conditions. 


Oxidation Through Coating Cracks 

Fig. 2 illustrates a typical oxidation void formed because of oxygen diffusion through a coating crack and 
subsequent reaction. A recent study by the authors outlines the microstructural features which distinguish this void 
from a processing void [4]. In general, an oxidation void is associated with a path for oxygen ingress (i.e., the crack 
in Fig. 2). At higher temperatures (>1000 °C) these oxidation processes are expected to be diffusion controlled. 
Diffusion control implies a lack of preferential oxidative attack and hence uniform attack of the carbon/carbon with 
a plane interface at the void boundary, as illustrated in Fig. 1. However, a close examination of the void in Fig. 2 
indicates some preferential attack along fiber edges and the resultant “pointed fibers,” which are characteristic of 
reaction control [5]. Thus the actual process is a mix of diffusion and reaction control. It will be shown that diffusion 
control dominates. 

There are numerous studies of carbon oxidation through cracks and fissures in the protective coating and 
associated models found in the literature. Medford [6] developed models of oxidation for RCC, assuming pathways 
for oxygen are the cracks and fissures in SiC. His model encompasses the following processes: 

(1) Diffusion of oxygen inward through the cracks and/or fissures in the SiC 

(2) Oxidation of SiC crack walls to form Si0 2 

(3) Oxidation of carbon/carbon — matrix, fibers, or both 

(4) Diffusion of CO outward through the cracks/fissures in the SiC 

Medford developed total mass change expressions based on these processes. The expressions were the diffusion 
equations with modifications for the area of reaction. First he examined the case below 982 °C, where passive 
oxidation of the SiC is thought to be negligible; then he examined the case above 982 °C, where passive oxidation of 
the SiC is thought to be significant. 

The general process of carbon oxidation in a nonreactive matrix or through a nonreactive coating has been 
explored by several investigators [7]-[ 12]. Carbon oxidizes to CO(g): 

C(s) + 14 0 2 (g) = CO(g) (1) 

However, thermochemically the equilibrium between 0 2 , CO, and C0 2 favors the products in the following reaction: 

1 /20 2 (g)+C0(g) = C0 2 (g) (2) 

Thus 0 2 (g) and CO(g) cannot coexist in similar amounts — they immediately will react to form C0 2 (g). However, 
the pairs of species in reaction (2), CO(g) with C0 2 (g) and 0 2 (g) with C0 2 (g), may coexist. Thus a two-step 
oxidation process is required to model carbon oxidation. The two reactions are 

1. At the carbon/gas interface: 
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C(s) + C0 2 (g) = 2CO(g) 


( 3 ) 


Note that C0 2 (g) is the oxidizing species. 

2. At a distance away from the carbon/gas interface: 

co(g) + y 2 o 2 (g) = co 2 (g) (4) 

In an earlier paper [12], the authors adapted this two-step diffusion-controlled oxidation model to cylindrical 
pinholes in the conversion coating. Actual oxidation damage was assessed with weight loss measurements and area 
measurement of oxidation voids below the pinholes, using cross sections. At 1400 and 1000 °C hemispherical voids 
were observed on a macroscopic scale, indicating uniform attack and diffusion control. The two-step model 
described the process reasonably well at 1400 and 1000 °C. However, at 600 °C only minimal oxidation is 
measured, and reaction control dominated. Hence the model did not fit at this lower temperature. 

In the present study, the two-step oxidation model is further extended to describe a half-cylinder trough which 
forms under a slot or craze crack in the SiC coating. First a series of experiments is performed using machined slots 
in the coating. This allows testing the model with well-defined oxygen pathway geometry. Then the study is 
extended to actual through-thickness craze cracks in the coating. Oxidation damage is assessed with cross-sectional 
views, interrupted or real-time weight loss measurements, and x-ray computed tomography scans. The kinetics of 
carbon consumption is then compared with the model. 


Experimental 

Characterization of Slots and Cracks in Test Specimens 

The tests described here were done on 1.91-cm-diameter RCC disks with a SiC coating on all sides. These disks 
were obtained from Lockheed-Martin Missiles and Fire Control in Dallas, Texas. Two sets of experiments were 
completed. Specimens with machined slots were used for the first set; specimens with naturally occurring craze 
cracks were used for the second set. Prior to oxidation, both types of specimens were characterized in order to 
quantify the area of exposed carbon. 

In the first set of specimens, artificial craze cracks with well-defined geometries were made with diamond 
blades (Keen Kut Products, Hayward, CA) of 0.25, 0.51, 0.76, and 1.02 mm thicknesses. These slots were cut to the 
SiC/carbon-carbon interface. These slotted specimens are shown in Fig. 3. One slot was cut in each specimen where 
the oxidation kinetics was followed with weight losses; however, two slots of different dimensions were cut in the 
specimens where the kinetics was followed using microstructure examination and/or x-ray computed tomography. 
Slot dimensions were carefully measured with a traveling microscope as well as directly from the cross sections and 
reported in the results section. Total coating thicknesses were measured from cross sections. Because of the 
variation in coating thickness, 20 random measurements were taken and averaged. 

The second set of specimens contained only the naturally occurring craze cracks. These were characterized with 
techniques similar to Yurgartis et al. [2]. Examination of the surface (plan view) did not reveal the complete crack 
pattern; however, grinding -300 pm of SiC clearly revealed the crack pattern, which was then traced, as illustrated 
in Figs. 4(a) and (b). Image analysis software (Adobe Photoshop CS2 (Adobe Systems, Inc., San Jose, CA) with 
Fovea Pro 4.0 (Reindeer Graphics, Asheville, NC)) was used to determine the total crack length per unit area. The 
view of a section, as shown in Fig. 5, allowed measurements of crack width at room temperature. Several such 
sections were examined, and an average of 10 measurements is reported. 

The measured values are listed in Table I and compared to the values reported by Yugartis et al. [2] and 
Buchanan and Little [3], “Specimen A” of Yurgartis et al. is SiC-coated carbon/carbon with no matrix inhibitor; 
“Specimen B” is SiC-coated carbon/carbon with a boron-based matrix inhibitor. It should be noted that the materials 
examined by Yugartis et al. [2] and Buchanan and Little [3] are different from the RCC materials in this study. Here, 
the RCC had a thicker SiC coating applied by a pack cementation process than the other specimens in which the 
coating was applied by a chemical vapor deposition process. Yugartis et al. examined cracks that do and do not 
transverse the entire coating thickness. In this study only through-thickness cracks are examined. Differences in the 
specific CTE of carbon/carbon and the applied coating would be expected to have a major effect on the crack 
parameters. Thus the differences listed in Table I are not unexpected. It is interesting that despite the differences in 
crack width and spacing, the ratio of the exposed area of carbon to the surface area of the sample is roughly the same 
for all materials. 
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Oxidation Treatments and Sample Analyses 

The samples were oxidized in air in either a box furnace or a vertical tube furnace. The kinetics of the box furnace 
exposures were followed with weight measurements every 0.5 h for a total of 2.5 h at 1200 °C. The craze crack 
specimens were suspended from a Cahn C-1000 (Cahn Instruments, Cerritos, CA) recording electrobalance by a 
sapphire fiber with a platinum wire basket into the vertical tube furnace. To avoid platinum silicide formation, the 
platinum wire was lined with alumina beads. Bottled air at 1 atm (101 325 Pa) was used at a flow rate of 1.67xl0 -3 
liter/s, and isothermal oxidation tests were performed at 1000, 1100, 1200, and 1300 °C. 

Selected samples were mounted and sectioned for optical microscopy to obtain measurements of oxidation void 
size. A specialized sample preparation procedure involved a premount of the entire sample in epoxy. A standard 
vacuum infiltration was done followed by a high-pressure application to fill the open porosity. Then the sample was 
sectioned, remounted in a standard metallographic mold, and polished with diamond pastes. This ensured clear 
microstructures and preserved the interface between the SiC and the carbon/carbon. 

X-ray computed tomography (CT) was used to provide additional images of the oxidation damage, without 
sectioning. This SmartScan Model 100 (CITA Systems, Inc., Pueblo, CA) system utilizes a Feinfocus FXE-160 
(COMET AG, Flamatt, Switzerland) microfocus x-ray source to produce very high-resolution imaging of samples, 
approaching 25 pm, in the CT mode of operation. The major hardware components of this system included the x-ray 
source, an area detector system, a five-axis object positioning subassembly, and a lead-lined radiation cabinet. A 
dual-processor computer system controlled the data acquisition and image processing. The slice plane thickness was 
120 pm per slice for these samples. Putting together slices electronically gave a three-dimensional view of oxidation 
damage. 


Oxidation Model Adapted to Slots and Cracks 

Fig. 6 is a schematic of the two-step diffusion-controlled model for oxidation through a pore or crack. Note that C0 2 
is the oxidizer and reaction (3) occurs at position x = 0 whereas reaction (4) occurs at position x = Xf. A list of 
symbols used in this report is in the appendix to aid the reader. 

The molar flux, J h of each species, z, is defined by 


J, = A eff 


r dc^ 

dx 


+ vr 


(5) 


The first term is the diffusive term and the second is the convective term. Here D, eff is the effective gas phase 
diffusivity, c, is the molar concentration of species i (i = CO, C0 2 , or 0 2 ), and x is the distance down the slot or 
crack. The average molar velocity vf ve is defined as 


5 >, Yj J i 



i 


Here Ct is the sum of concentrations. Equation (5) becomes 


Ji=Df 


r dc A 






( 6 ) 


(7) 


The boundary conditions are shown in Fig. 6(b): 


At x = 0 : c co = c£ 0 c COl = c° 02 
At x = Xf : c c q 2 = cq 2 = cco = 0 


( 8 ) 
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At x = L\ c 0l = c £ 2 c C o 2 = 0 


As shown in Fig. 6(b) L is the total crack depth, which is the coating thickness. 

Consider the fluxes in Region 1 (Fig. 6(a)). Flere = -Jqq and hence the convective term in equation (7) is 

zero. The fluxes of 0 2 (g) and C0 2 (g) in this region are thus calculated as 


J l o 2 =-°o 2 


dcp 2 

dx 


- °<h c o 2 

l L ~ x f ) 


(9a) 


co- 


= D, 


co 2 


dc. 


co 2 


dx 


D C0 2 c co 2 
( L ~ x f) 


(9b) 


A determination of the C0 2 flux emerging from the crack or slot, Jqq , is the ultimate goal, because it is a measure 

of carbon consumption. First the quantity Xj i s needed. This can only be obtained as the ratio Xf/L, which is derived 
from the following equations (eqs. (10) to (17)). The two fluxes in equation (9) are equal and opposite. From this 
fact, one can deduce the concentration gradients in Fig. 6(b). The concentration of C0 2 at x = X/is 


c 


co 2 


D Oi c o 2 
d co 2 


(10) 


Equation (7) for CO and C0 2 in Region II (Fig. 6(a)) is more complex, as now there is a convection term since the 
stoichiometry of reaction (3) requires 


,/H - -2J 11 

,y co ZJ co 2 


Thus equation (7) for C0 2 becomes 


J" =~D, 


co 


co 2 


dc , 


co 2 


dx 


c CO^J C( 


2 l/ C0 2 


Ct 


Rearranging to separate variables gives 


^'o 2 ^ 


co- 


C T 


dc, 


co 2 


c T +c co , j 

This can be integrated from c C q 2 = at x = xy and c co , = c® 0 at x = 0 : 




ct +c. 


co 2 


ct +c. 


co 2 


This can be simplified further using equation (10) and the fact that c£ 0 » c° Q 


(11) 


( 12 ) 


(13) 


(14) 


J 


II 

co 2 


-Dqo 2 CT 

x f 


In 


D, 


1 + - 


°2 C Q 2 

Dqo 2 c t 


(15) 


6 



A mass balance between CO and O 2 gives 


jii =-2 
J CO o. 


(16) 


Combining this with equation (11), 


J 


11 

C0 2 



When equation (15) is combined with equations (17) and (9a), 


Solving for i/x/- gives 


j 11 

J co 2 


°o 2 c b 2 

[ L ~ x f) 


-A 


'co 2 c t 


In 


x f 


1 + - 


^°2 C 0 2 


A 


'co 2 c t 


L { ( °co 2 c ( j 2 

In 

fi+^M 

Xf D Ql c T 


Dco 2 c t J 


(17) 


(18) 


(19) 


Next the gas phase diffusivities, A.n, > of 0 2 in N, and C0 2 in N 2 must be estimated. This is done with the 

Chapman-Enskog correlation [13], which yields diffusivity in cm 2 /s when the variables in the expression below are 
expressed in the appropriate units: 


A n , =0.0018583 


1 

M 7 




N t/2 f 

fill 


2 J 


V 


Ho 2 , Q 

;,N, J 


(20) 


Here M, is the molecular weight (g/mol) of species i, T is the absolute temperature (K), P is the pressure (atm), 
o, n 2 * s the average of the molecular diameter (A) of species i and N,, and Q;n 2 is the collision integral (unitless). 

The later two were determined from the tables of Svehla [14] and Sherwood et al. [15]. Table II lists some 
calculated diffusivities. The value of Ax/is also shown in Table II and is constant with temperature. 

For narrow cracks, the contribution of Knudsen diffusion must also be considered. In the Knudsen diffusion 
regime, molecule-wall collisions dominate over molecule-molecule collisions. The effective diffusivity, D e " is a 
combination of the ordinary diffusivity, D ord , and Knudsen diffusivity, D K . The Knudsen diffusivity can be estimated 
for a pore as [16] 




(21) 


Here q is the pore radius in centimeters. Equation (21) is for a pore, not a crack, but it is still suitable for this general 
discussion of Knudsen effects. The combined effects of ordinary and Knudsen diffusion can be written as [17] 


1 


1 1 
- + - 


Z)eff £)ord [)k 


(22) 


When the mean free path is greater than 10 crack widths, then Knudsen diffusion is the dominating flow mechanism 
[18]. For the 1 atm total pressure used in this study, the mean free path is -3.3x10 m. The machined slots are much 
larger than this, and even the craze cracks (-1 ,28x 1 0 m) are larger than the mean free path. This becomes apparent 
by comparing D° rd /D c{t to the pore diameter, as shown in Fig. 7. The level parts of the plots show A ff = A >ld . At the 


7 



1 atm (101 325 Pa) pressure used in this study, the average crack width indicates a very small Knudsen contribution; 
however, at lower pressures which are closer to actual reentry pressure, the Knudsen contribution must be 
considered. 

Now with values for Xf/L and the gas phase diffusivities carbon consumption can be calculated. As noted, this is 
simply the molar flux of COdg) emerging from the crack, as given by equations (9b) and (9a): 


J 


I 

co 2 


d co 2 c C0 2 

l - L ~ x f ) 


D 0i c 0 2 


1- 




(23) 


The denominator is written in the above form, since Xf/L is now a known quantity. The goal is to estimate the radius 
of the cavity formed below a crack, as shown in Fig. 6. 

The flux of C0 2 (g) (eq. (11)) can be converted to a mass loss of carbon as 


d W c 
d t 


= M cJ ho 2 wl = 


McDo 2 Po 2 wl 


RTL 


1- 




(24) 


Here W c is the weight loss of carbon, M c is the molecular weight of carbon (12 gm/mol), R is the gas constant, w is 
the crack width, and / is the crack length. Note that the wl term would be omitted for weight loss per unit area. 

The weight loss with time can be related to the volume change with time of the cavity at the base of the crack. 
Assume a rectangular crack with a half-cylinder cavity of radius r growing at its base. In practice the nonideality of 
the crack is accounted for with a “tortuosity factor” [19]. For this approximation a tortuosity factor of 1 is assumed. 
The weight loss can be related to the radius of the growing half-cylinder cavity as 


dW c dV 
— — = P — 
d 1 dt 



(25) 


Here r is the radius of the half cylinder and p is the density of carbon/carbon taken to be 1.362 g/cm 3 [20]. Equating 
equations (24) and (25) and separating variables gives 


nr 2 l 


M c Do 2 P 02 wldt 


pRTL 


1- 


YY 

k l J 


(26) 


Integrating and solving for r results in the following: 


2M c Dq 2 Pq 2 wt 



( ,, Y 


npRTL 

H 

x f 



\ L 



(27) 


Note the crack length cancels out and does not appear in the final expression. Equations (24) and (27) are the key 
equations used to calculate oxidation damage in subsequent sections. 

An important consideration is the variation of the crack width with temperature. This is minor for the machined 
slots, but significant for the craze cracks. The actual coating is not dense SiC, but rather converted carbon-carbon 
fibers and carbon matrix material. For this approximation the accepted coefficient of thermal expansion (CTE) for 
SiC as ~6. 1 x 1 0 fl K _1 is used [21]. The thermal expansion along a carbon fiber is taken as 1 x 1 0 6 K _l [22]. This puts 



the coating in tension on cooldown, which leads to the cracks during processing. However, it leads to compression 
and possible crack closing during heatup to temperatures above the processing temperature. The extent of crack 
closure can be estimated as 


W RT = Aas(T P - T) (28) 

Here W Rr is the crack width at room temperature, Aa is the difference between the CTE of SiC and that of the 
carbon/carbon, s is the crack spacing, T P is the processing temperature (taken as -1650 °C), and T is the operating 
temperature. Since the crack width at room temperature is 12.8+1 .4 pm, the crack spacing is calculated to be 1.5±0.2 
mm from the equation above, reasonably close to that measured given the approximation of the CTEs. This 
calculated crack spacing is used in the above equation to estimate crack width at the oxidation temperature. At 1200 
°C the crack width was calculated to be 7.1 pm. 

Oxidation along the SiC crack wall may affect the rate of carbon/carbon oxidation. As will be discussed, the 
samples gain considerable weight due to passive oxidation. Microscopy indicates that thin films of Si0 2 form on 
both the internal and external surfaces of the SiC coating. Fig. 8 shows the mouth of a crack and only a thin (-0.25- 
pm) film of Si0 2 . 

On the basis of the thin films of oxide observed, the oxidation rate constants for high-purity chemically vapor 
deposited (CVD) SiC [23] are used. The results shown in Table III are somewhat lower than the thicknesses of the 
oxide film on both the slots and cracks. It is probable the lower purity SiC in the conversion coating has an oxidation 
rate constant of an order of magnitude or so greater than that of CVD SiC. Such behavior is well documented [24]. 
Nonetheless images such as those in Fig. 8 indicate that under the conditions of this study, crack closure due to wall 
oxidation can be neglected. 

Results and Discussion: Measured Oxidation Kinetics and Comparison With Model 

Oxidation Below Machined Slots in SiC Coating 

Initially, two specimens were oxidized at 1200 °C and examined with both CT and optical microscopy. One specimen 
had nominal 0.25- and 0.51-mm slots machined in the face. The other specimen had nominal 0.76- and 1.01-mm slots 
machined in the face. Cross-sectional micrographs for the cavities formed below the nominal 0.25-mm slots are shown 
in Figs. 9(a) and (b). Oxidation damage below all slots exhibited this type of half-cylinder morphology. Figs. 10(a) and 
(b) show the x-ray CT information. 

The data in Figs. 9 and 10 make several important points. The similarity between the CT images and the actual 
optical microscopy images give credibility to the use of CT imaging to precisely size and characterize the three- 
dimensional morphology of this type of oxidation damage [25]. Both the optical microscope images and the CT 
images indicated the oxidation damage occurs in the form of a half cylinder. On a macroscopic scale this indicates 
even attack in all directions, supporting diffusion control of the oxidation process and the use of a half cylinder in 
the model. 

Using image analysis software the areas of the oxidation cavities were measured. Approximating these areas as 
semicircles allows calculation of effective radii. The experimental results are compared to the predictions from 
equation (27) in Table IV. The error in the predicted oxidation slot radius values is due to uncertainty in the 
measurement of the slot widths. Agreement between the measured and predicted cavity radius is within a factor of 
0.45 to 0.62. All the measurements are consistently lower than the predictions, which is likely due to carbon 
consumption beyond the apparent boundaries of the half cylinder, as seen clearly in Fig. 9(b). 

A second set of experiments was done to follow the kinetics of oxidation below machined slots with weight loss 
as the indicator of oxidation damage. In this case specimens with only a single slot were used. Kinetics were 
followed using interrupted weight loss measurements. A typical kinetic curve is shown in Fig. 11. Note there is an 
initial weight gain, which is likely passive oxidation of the SiC, followed by a weight loss. The passive oxidation of 
the SiC will be discussed in the following section. The rates of weight loss are the carbon oxidation rates per the 
exposed area of carbon, which was determined from the slot dimensions. Table V compares the measured rates with 
those predicted. Rates are normalized per unit area of exposed carbon, and hence the predictions are all the same for 
the given temperature of 1200 °C. The error in the predicted oxidation rate values is due to uncertainty in the 
measurement of slot depth (coating thickness). The error in the measured rate values is due to uncertainty in the 
measurements of slot width and length. Note that these errors make the measured rates for three of the four smaller 
slots within the limits of the predictions (Table V). For the larger width slots (-1 mm), the measured rates are -0.33 
of the predicted. The reasons for this are not clear. 
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Oxidation Below Craze Cracks in the SiC Coating 

Having established the validity of the model with the machined slots, it is now appropriate to proceed to the more 
realistic situation of craze cracks in the SiC coating. Now there is a less well-defined geometry, as illustrated with 
the crack pattern in Fig. 4 and crack width measurement in Fig. 5. Nonetheless these measurements can be used to 
estimate approximate values for the diffusion model. 

Oxidation kinetics for these specimens were followed with a recording thermobalance that enabled continuous 
weight change measurements. Results for 1000, 1100, 1200, and 1300 °C are given in Figs. 12(a) to (d). These data 
are reported in weight change per geometrical surface area. As will be discussed, the geometric surface area must be 
modified for actual reactive area of oxidation, which is different for each of the two processes occurring. 

The oxidation kinetics suggest a combined parabolic/linear rate law. The initial weight gain observed at the 
three higher temperatures is likely due to the passive oxidation of SiC, which is parabolic: 



(29) 


Here 



is the specific weight change (weight per unit area), k p is the parabolic rate constant, and t is the time. 


The specific weight change over the surface area A of the SiC is taken as the geometric area of the sample disk. The 
weight loss after ~0.5 h is linear with time and most likely due to carbon oxidation through the craze cracks. It is 
described by 



(30) 


Here k/ is the linear rate constant. Note that A' is the exposed area of carbon/carbon. This is the crack width at 
temperature multiplied by the total crack length for a particular sample. 

The changes in the kinetic curves with temperature support the interpretation of the process as passive SiC 
oxidation coupled with C/C oxidation. At 1000 °C the initial parabolic regime was not observed; however, at 1200 
and 1300 °C it was quite evident. This is due to the fact that passive oxidation of SiC shows a strong dependence on 
temperature [23]. Note, however, that the linear rate constant is relatively temperature independent. Diffusion- 
controlled oxidation of carbon shows only a weak temperature dependence, a combination of the T ] 2 dependence of 
the diffusion coefficient and the T l dependence of the crack closure due to thermal expansion. This dependence may 
not be evident in the limited data set in this study. 

The combined passive oxidation of SiC/linear oxidation of carbon rate law is written as 



(31) 


is the measured weight change per unit of geometrical area, k p is the parabolic rate constant, t is 

is the linear rate constant. This expression is similar to that used to describe the combined 
oxidation/vaporization of chromium alloys [26] and SiC [27], but this is somewhat simpler in that there are two 
independent processes occurring. 

The dashed lines in Figs. 12(a) to (d) are fits of equation (31) to each kinetic curve. From these fits the parabolic 
and linear rate constants can be derived. Fits to the linear region are good; fits to the parabolic region are only 
approximate because of such factors as changing surface areas of passive oxidation and passive oxidation internal to 
the coating, which are not included in equation (31). The parabolic and linear rate constants extracted from the fit to 
the kinetic data are given in Table VI. The parabolic rate constant was taken directly from the fit. Note that the listed 
parabolic rate constants are many orders of magnitude greater than those measured for high-purity SiC [23]. This is 
due to extensive internal oxidation of the porous SiC conversion coating. Measurements of the SiC conversion 
coating density have been reported to be 2.114 g/cm 3 [20]. When compared with the theoretical density of 3.217 
g/cm 3 for SiC, this indicates a high amount of porosity. This porosity leads to the large extent of internal oxidation. 



the time, and ki 
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The linear rate constant shown in Table VI was taken from the fit and multiplied by the geometrical surface 
area. This gives a rate of carbon oxidation per area of exposed carbon/carbon and allows direct comparison to the 
calculated rate in Table VI. The predicted rates were obtained from equation (24). 

The cracks in the SiC form a complex geometry. The oxidation cavities are an irregular shape as seen in the 
micrograph in Fig. 2. Fig. 13(a) shows the location of two CT scans, and Fig. 13(b) shows the results of the CT 
scans illustrating the varying size and shape of the oxidation cavity below a craze crack. The actual oxidation 
cavities are darker and are outlined in white in Fig. 13(b). 

Despite this complex geometry of both the craze crack and the resultant oxidation damage, average values of 
crack width and length can be put into the oxidation model to estimate the rates. These are shown in Table VI. 
Experimental rates were found to be both greater and less than the predicted rates for the temperature range of 1000 
to 1300 °C. As in the case of the slots, the error associated with the predicted oxidation rate values is due to 
variation in the measurement of crack depth (coating thickness), and error associated with the measured rate values 
is primarily due to uncertainty in the crack length and width measurements. The predicted and measured rates show 
agreement within experimental error, which is remarkable given the irregular geometry of the cracks and 
corresponding oxidation damage as well as the approximations used in this study. 


Conclusions 

The kinetics of subsurface oxidation of SiC-protected reinforced carbon/carbon has been studied in the diffusion 
control regime (1000 to 1300 °C) in air. Specimens with machined slots in the SiC coating as well as specimens 
with natural occurring craze cracks in the SiC coating were tested. The geometry of these oxidation pathways was 
characterized with the aim of obtaining the exposed area of the carbon/carbon. For the craze cracks, crack closure 
due to shrinkage was considered. Oxidation of the SiC walls was shown to be negligible under the conditions of this 
study. Oxidation damage of the carbon/carbon was assessed with weight changes, CT scans, and optical microscopy 
of sections. The CT scans confirm the idealized half-cylinder shape of the oxidation cavity in the case of a slot. A 
two-step diffusion control model for carbon oxidation is adapted to describe the growth kinetics of this half cylinder. 
Oxidation rates were first measured for the well-defined geometry of the slotted specimens. In the case of optical 
measurements of oxidation, the error associated with the predicted values of oxidation rate were primarily due to 
uncertainty in the measurements of the slot width, and the error in the measured rate values were due to uncertainty 
in measuring the size of the oxidation cavities. For these optical measurements, measured rates were consistently 
lower than the predicted by a factor ranging from 0.45 to 0.62. For the weight loss measurements of oxidation, errors 
in the predictions were primarily due to errors in measuring the slot depth (coating thickness) and errors in the 
measurements were due to errors in measuring the slot width and length. Most of the predictions were within 
experimental error of the measurements, except for the largest slots (~1 mm width), where the measurements differ 
by a factor of -0.33 from the predictions. Proceeding next to the actual craze cracks, the oxidation cavity was found 
to be highly irregular and approximated by a half-cylinder shape. Weight loss was used to measure oxidation, and 
again, errors in the predictions were primarily due to errors in measuring the slot depth (coating thickness) and 
errors in the measurements were due to errors in measuring the crack width and length. For oxidation below craze 
cracks, the predictions are generally within the experimental error of the measurements. 
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Appendix — Symbols 


A 

A’ 

c 

D 

J 

k 

L 

l 

M 

P 

q 

R 

r 

s 

T 

t 

v 

W 

w 

X 

x f 

a 

Q 

a 


area of SiC for oxidation 

area of exposed carbon-carbon for oxidation 

molar concentration 

diffusivity 

molar flux 

rate constant (weight 2 /(unit area) 2 -time) 
depth of crack or coating thickness 
crack length 

molecular weight (g/mol) 
total pressure in atmospheres 
pore radius 
gas constant 

radius of half-cylinder cavity 
crack spacing 
temperature 
time 

molar velocity 
weight loss 
crack width 

depth of crack coordinate with origin at SiC/C-C interface 
depth to C0/C0 2 (region II) and C0 2 /0 2 (region I) boundary 
average collision diameter between species i and N 2 (angstroms) 
collision integral 
coefficient of thermal expansion 


Subscripts: 

/ 


/ 

P 

P 

T 


indicates boundary of regions I and II 

species 

linear 

processing 

parabolic 

total 


Superscripts: 


ave 

average 

eff 

effective 

K 

Knudsen 

L 

indicates position L 

ord 

ordinary 

0 

indicates position 0 

I 

region I 

11 

region 11 

* 

indicates position x f 
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TABLE I Measured crack parameters for SiC-coated carbon/carbon materials 


SiC-coated carbon/carbon 
material 

Coating 

thickness, 

mm 

Crack length per 
unit area, 
mm 1 

Crack spacing, 
mm 

Crack width, 
pm 

Crack area per 
unit area 

Present study, RCC 

0.78+0.14 

0.33+0.04 

3. 4+0. 9 

12.8+1.41 

4.2x1 0“ 3 

Yurgartis et al. [2], Specimen A 

0.3024+0.03866 

2.985 

0.6247+0.2233“ 

1.72+1.41 

5.1xl0“ 3 

Yurgartis et al. [2], Specimen B 

0.2788+0.02853 

3.555 

0.6517+0.3356“ 

2.02+1.47 

7.3xl0“ 3 

Buchanan and Little [3] 

— 

— 

.... 

— 

Q-SlxlO" 4 


a Through coating thickness. 
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TABLE II Diffusivities of CLIg) and CCLfg) in N 2 (g) and 
position of boundary between regions I and II (Fig. 6(a)). 
[Total pressure is 1 atm (101 325 Pa).) 


Temperature, 

°C 

0 2 ,N 2 

cm 2 /s 

Ax> 2 ,n 2 ? 
cm 2 /s 

Ltxf 

1000 

2.38 

1.86 

2.13 

1100 

2.70 

2.11 

2.13 

1200 

3.03 

2.38 

2.13 

1300 

3.38 

2.65 

2.13 
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TABLE III Estimated oxide film thicknesses on slot or crack walls of SiC coating 


Temperature, 

Rate constant 

Oxide thickness 

°C 

1 atm 0 2 , a 

0.21 atm 0 2 , 

0.5 h in 0.21 atm 0 2 , 

2.5 h in 0.21 atm 0 2 , 


pm 2 /h 

Hm 2 /h 

pm 

pm 

1200 

1.82xl0“ 2 

3.78xl0' 3 

4.35xl0~ 2 

9.72xl0“ 2 

1300 

3.76xl0 -2 

7.90x1 O' 3 

6.28xl0' 2 

0.14 


“From reference 23. 
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TABLE IV Measured and calculated radii of oxidation 
cavity for SiC-coated RCC specimens with machined slots 
[In static air at 1200 °C.] 


Width of slot, 

Measured 

Calculated radius, 3 

mm 

radius, 

r. 


mm 

mm 

0.5 h exposure 

0.26+0.02 

0.99+0.07 

1.82+0.07 

0.53+0.02 

1.22+0.03 

2.60+0.05 

0.81+0.02 

1.51+0.03 

3.32+0.04 

1.08+0.01 

1.51+0.03 

3.77+0.02 

2.5 h exposure 

0.44+0.08 

2.38+0.01 

3.80+0.30 

0.59+0.08 

2.60+0.16 

4.40+0.26 

1.05+0.07 

2.96+0.04 

5.87+0.20 


a From equation (27). 
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TABLE V Oxidation rates for SiC-coated 
RCC specimens with machined slots 
[At T= 1200 °C in air.] 


Slot width, 
w, 
mm 

Slot length, 
/, 

mm 

Weight loss rate, 
mg/h-mm 2 

Measured Calculated 11 

0.484+0.062 

8.177+0.923 

34+6 

55+10 

0.312+0.033 

8.284+0.63 

48+6 

55+10 

0.560+0.013 

7.094+0.132 

48+1 

55+10 

0.466+0.04 

7.114+0.019 

57+5 

55+10 

1.116+0.024 

7.602+0.101 

18+1 

55+10 

0.963+0.0041 

7.624+0.039 

21+1 

55+10 


a From equation (24). 
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TABLE VI Measured and calculated rates of RCC oxidation through craze cracks 


Sample temperature, 
°C 

Geometrical 
surface area, 
mm 2 

Total crack 
length, 3 
mm 

Crack width, 
pm 

Area of carbon 
exposed by 
craze cracks, b 
mm 2 

Calculated 
linear rate, c 
mg/mm 2 -hr 

Measured 
linear rate, c ’ 
mg/mm 2 -hr 

Measured 
parabolic rate, e 
mg 2 /mm 4 -hr 

Room temperature 
1000 

850.2+10 

281+34 

12.8+1.4 

10.2+2.7 

2.87+0.8 

50+9 

31+9 


1100 

935.5+10 

309+37 

8. 7+2. 3 

2.67+0.8 

53+9 

31+9 

3.94x10^ 

1200 

946.8+10 

312+37 

7.1+1. 8 

2.21+0.6 

55+10 

53+15 

2.79x1 0- 2 

1300 

940.2+10 

310+37 

5.5+1 .4 

1.71+0.5 

57+10 

67+20 

3.25xl0~ 2 


“Equal to (geometrical area x crack length)/unit area. 

b Equal to (total crack length) x (crack width). 

c From equation (24). 

d Area is exposed area of carbon/carbon. 

c Area is geometrical surface area. 
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Fig. 1. — Reinforced carbon/carbon (RCC) sample. 

Fig. 2. — Optical micrograph showing oxidation void in RCC formed below crack in SiC coating at 1100 °C and 0.0066 atm (667 
Pa) for 1 h. 

Fig. 3. — Machined slots in RCC buttons. 

Fig. 4. — Crack pattern on RCC disk, (a) Optical micrograph, (b) “Skeleton” trace of cracks. 

Fig. 5. — Crack in SiC conversion coating. 

Fig. 6. — Two-step diffusion control model [12], (a) Schematic of fluxes and boundaries, (b) Relative concentrations of gases. 

Fig. 7. — Plot to show relative contribution of Knudsen diffusivity to effective diffusivity. The vertical line shows average crack 
width. 

Fig. 8. — Oxide film on surface of SiC coating (0.5 h in air at 1200 °C) at crack opening, (a) Electron micrograph, (b) Oxygen 
map. 

Fig. 9. — Oxidation below nominal 0.25-mm slots in SiC coating, (a) 0.28-mm slot in air at 1200 °C for 0.5 h. (b) 0.37-mm slot in 
air at 1200 °C for 2.5 h. 

Fig. 10. — Computed tomography (CT) of machine-slotted SiC-coated RCC oxidized for 0.5 h in air at 1200 °C. (a) CT sample 
locations, (b) Three CT slices. Slot near notch on left is 0.76 mm, and slot on right is 1.03 mm. (c) Three-dimensional 
reconstruction of disk. 

Fig. 11. — Oxidation kinetics determined from interrupted weight loss measurements for RCC specimen with 0.484-mm slot 
width at 1200 °C in air. 

Fig. 12.— Oxidation kinetics for RCC in air. (a) 1000 °C. (b) 1 100 °C. (c) 1200 °C. (d) 1300 °C. 

Fig. 13. — Computed tomography (CT) of craze crack of SiC-coated RCC oxidized for 2.5 h in air at 1 143 °C. (a) Location of CT 
slices, (b) Two CT slices. 
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